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Abstract. 

Aims: In this work the luminosity function of QSOs is measured in the redshift range 3.5 < z < 5.2 for the absolute magni- 
tude interval -21 < M 145 < -28. Determining the faint-end of the luminosity function at these redshifts provides important 
constraints on models of the joint evolution of galaxies and AGNs. 

Methods: We have defined suitable criteria to select faint QSOs in the GOODS fields, checking their effectiveness and com- 
pleteness in detail. A spectroscopic follow-up of the resulting QSO candidates was carried out. The confirmed sample of faint 
QSOs is compared with a brighter one derived from the SDSS. We used a Monte-Carlo technique to estimate the properties of 
the luminosity function, checking various parameterizations for its shape and evolution. 

Results: Models based on pure density evolution show better agreement with observation than do models based on pure lumi- 
nosity evolution. However, a different break magnitude with respect to z ~ 2.1 is required at 3.5 < z < 5.2. Models with a 
steeper faint-end score a higher probability. We do not find any evidence for a bright-end flattening at redshift z > 3.5. 
Conclusions: The estimated space density evolution of QSOs indicates a suppression of the formation and/or feeding of su- 
permassive black holes at these redshifts. The QSO contribution to the UV background is insufficient for ionizing the IGM at 
3.5 < z < 5.2. 

Key words, quasars: general - galaxies: active - cosmology: observations 



& ' 1 . Introduction the dark matter halo (DMH) hierarchical sequence for the col- 

lapse of baryonic mat ter, a phenomenon also known by the 



Once considered peculiar and exotic objects, quasars (QSOs) term downsizing (e . g. khanate et alJl2004HCroton et all 



have been recognized in recent years as an important and | B ower et all EoOft jFontanot et alJ bOOft iMenci et al. I 



2006; 



- - - - - - - - - 2006; 

possibly a necessary phase of galactic evolution, at least for Scanna meco & Ohll2004 . The depth of the potential wells in 

spheroids (Danese et al. 2003). In particular, models of a joint different DMH i s a key factor in determining the timescale for 

evolution of galaxies and QSOs (e.g.| Granatoetalj200lH2004 star formation, accretion, and eventually for the effects of feed- 

| Hopkinsetalj2005H2006|) reproduce the properties of present- back ^ setting the pace for AGN activky mroughout cosmic his . 

day elliptical galaxies assuming that their vigorous star forma- tor y 
tion at high redshift is quenched by the feedback (i.e. galactic 

winds). This process is triggered by the QSO activity fed by the Great progress has recently been achieved thanks to the 
accretion of matter onto a supermassive black hole (SMBH) at large amount of data coming from major observational pro- 
file center of the galaxy itself (see also Monaco & Fontanot grams such as the Two Degree Field QSO Redshift Survey 
2005). In this way, the properties of high redshift QSOs and, (2QZ, Croom et al. 2004, www.2dFquasar.org) and the Sloan 
in particular, their luminosity function (LF) are fundamental Digital Sky Survey, with the third edition of the Quasar Catalog 
in understanding the phenomena driving galaxy formation. It (DR3QSO, Schneider et al. 2005). The 2dF-QRS contains 
has been hypothesized that the feedback from AGN inverts m0 re than 23000 QSOs typically at redshift z < 2.1, while the 

DR3QSO lists 46420 QSOs up to z ~ 5.4. However, at high 

Send offprint requests to: F. Fontanot redshift the SDSS is sensitive only to the most luminous QSOs 
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(M 145 £ -26.5 1 ). The faint-end of the high-z QSO LF, which 
plays a key role in comparing different predictions of the for- 
mation and evolution of galaxies, is left almost unconstrained 
due to the challenging magnitude limit. 

The recent Great Observatories Origins Deep Survey 
(GOODS) provides an important database for studying QSOs 
with M145 ^ —21, thanks to the depth of its optical observa- 
tions and to its multi-wavelength nature. GOODS is a deep 
survey, not a wide one, but it is much larger than most pre- 
vious deep HST/WFPC2 surveys, covering 320 arcmin 2 , 32 
times the combined solid angles of the Hubble Deep Field- 
North and South, and four times larger than their combined 
flanking fields. The GOODS-ACS optical data were analyzed 
in a previous paper (Cristiani et al. 2004 hereafter paper I) in or- 
der to select reliable QSO candidates and make a first estimate 
of their space density in the redshift interval 3.5 < z < 5.2. The 
selection was carried out by defining suitable optical selection 
criteria based on magnitude limits and color criteria, and then 
by matching the optical candidates with Chandra X-ray surveys 
(Alexander et al. 2003; Giacconi et al. 2002). 

In this paper we take advantage of the SDSS 3 rd Data 
Release, of a new analysis and of a more complete spectro- 
scopic follow up of the GOODS data in order to determine 
the QSO space density at 3.5 < z < 5.2, down to the faint- 
end of the LF (M145 < -21). As in paper I we adopt a def- 
inition of QSOs encompassing all objects with strong, high- 
ionization emission lines and M145 < -21, including conven- 
tional, broad-lined (type-1) QSOs, and narrow-lined, obscured 
(type-2) QSOs. Throughout the paper we adopt a flat universe 
with cosmological parameters Q m = 0.3, Oa = 0.7, and Hubble 
constant Hq = 70Km s~ l Mpc~ l . Magnitudes are in the AB sys- 
tem. 

In Sect. 2 we present the databases, the selection of the 
QSO candidates, and the spectroscopic follow up; in Sect. 3 
we present our algorithm for the determination of the LF; in 
Sect. 4 we discuss the results in terms of implications for cur- 
rent models of galaxy formation. 

2. The database 

2.1. Bright QSOs from the SDSS 

The bright quasar sample we use in this work was extracted 
from the DR3QSO. The main SDSS catalogue covers an area of 
about 5282 deg 2 in photometry and 4188 deg 2 in spectroscopy. 
Photometric catalogues were compiled with observations in u, 
g, r, i, z bands (Fukugita et al. 1996). Several selection crite- 
ria were tailored in order to select QSOs candidates in SDSS 
at different redshifts (see e.g. Fan et al. 2001). The selection 
criteria have been changing significantly over time (Schneider 
et al. 2005). Richards et al. (2002, hereafter R02) present the 
updated version of these selection criteria applied to select suit- 
able QSO candidates for spectroscopy. We refer to these crite- 
ria in the following sections (namely to eq. 6 and 7 in R02). In 
particular we are interested in the selection of high-z objects. 
Due to the strategy of the SDSS observations, the DR3QSO is 
a sample with an incomplete follow-up. 

1 Absolute AB magnitude at a wavelength of 145 nanometers 



Fan et al. (2003) have compiled a complete sample of high- 
z QSOs out of the SDSS commissioning data; however, we de- 
cided to use DR3QSO, in order to have a larger sample of ob- 
jects selected in a larger area of the sky. The price of this choice 
consists our being compelled to make a more complex analy- 
sis to take the incompleteness of the DR3QSO follow-up into 
account for our purposes (see Sect. 3). There are 656 QSOs 
in DR3QSO that have a redshift between 3.5 < z < 5.2 and 
that satisfy R02 criteria. We considered those objects in the 
SDSS photometric catalogue that satisfy the R02 criteria and 
we compared them to the corresponding SDSS spectroscopic 
catalogue. This way, we have estimated that only ~ 21% of the 
candidates satisfying the R02 criteria have a spectrum. 

2.2. Faint QSOs from the GOODS 

The GOODS covers an area of 320 arcmin 2 , subdivided into 
two 160 arcmin 2 sub-fields centered on the Chandra Deep 
Field-South (CDF-S) and Hubble Deep Field-North (HDF- 
N). The optical data (B445, V660, '775, Z&50 bands) were ob- 
tained with the Advanced Camera for Surveys (ACS) onboard 
HST in the framework of the GOODS/ACS survey described 
in Giavalisco et al. (2004a). The catalogues used in this pa- 
per, prepared using the SExtractor package (Berlin & Arnouts 
1996), are based on the version vl.O of the reduced, calibrated, 
stacked, and mosaiced images acquired with HST and ACS as 
part of the GOODS ACS Treasury program 2 . The catalogues 
are z-band-based, that is, source detection was made using the 
z-band images. The magnitude limits in the four bands are 
27.50 (B445 band), 27.25 (V 66 o band), 27.00 (/ 77 5 band), and 
26.5 (Z850 band) at S/N 10 for point sources. 

The HDF-N and CDF-S fields were observed in the X- 
rays with Chandra for 2 Ms and 1 Ms, respectively (Alexander 
et al. 2003, hereafter A03; Giacconi et al. 2002, hereafter 
G02), providing the deepest views of the Universe in the 0.5- 
8.0 keV band. The X-ray completeness limits over m 90% 
of the area of the GOODS fields are similar, with flux lim- 
its (S/N= 5) of * 1.7 x 10~ 16 erg cirT 2 s" 1 (0.5-2.0 keV) 
and * 1.2 x 10~ 15 erg cirT 2 s _1 (2-8 keV) in the HDF- 
N field and * 2.2 x 10~ 16 erg cirT 2 s" 1 (0.5-2.0 keV) and 
* 1 .5 x 10~ 15 erg cirT 2 s 1 (2-8 keV) in the CDF-S field (A03). 
The sensitivity at the aim point is about 2 and 4 times better for 
the CDF-S and HDF-N, respectively. As an example, assuming 
an X-ray spectral slope of F = 2.0, a source detected with a 
flux of 1.0 x 10~ 16 erg cirT 2 s _1 would have both observed and 
rest-frame luminosities of 8 x 10 42 erg s _1 and 3 x 10 43 erg s _1 
at z — 3 and z = 5, respectively (assuming no Galactic absorp- 
tion). A03 produced point-source catalogues for the HDF-N 
and CDF-S and G02 for the CDF-S. 

2.2.1 . The color selection 

The selection of the QSO candidates was carried out in the 
magnitude interval 22.25 < Z&50 < 25.25. To avoid the lower 
quality zones at the borders of the GOODS ACS mosaics, we 



2 see, for details http://www.stsci.edu/science/goods/ 
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carried out the selection in a slighty reduced area of 157.1 
arcmin 2 in the HDF-N and 156.4 arcmin 2 in the CDF-S. 

Expected QSO colors in the ACS bands were estimated 
as a function of redshift using a template of Cristiani & Vio 
(1990, hereafter CV90) for the QSO spectral energy distribu- 
tion (SED) convolved with a model of the intergalactic medium 
(IGM) absorption (see Appendix). The same scheme will be 
adopted in the following for the estimate of expected colors in 
the different SDSS passbands, unless otherwise specified. 

In paper I, four optical criteria had been tailored to select 
QSOs at progressively higher redshift in the interval 3.5 H z ^ 
5.2.: 

i-z < 0.35 AND 1.1 <B- V< 3.0 AND V -i< 1.0 (1) 
i-z < 0.35 AND B - V > 3.0 (2) 
i-z < 0.5 AND B - V > 2.0 AND V-i> 0.8 (3) 
i-z < 1.0 AND V-i> 1.9. (4) 

As already shown in paper I, these color selection criteria are 
expected to be most complete and reliable at z > 4 (correspond- 
ing to criteria 2-A). These criteria select a broad range of high- 
z AGN, not limited to broad-lined (type-1) QSOs, and are less 
stringent than those typically used to identify high-z galaxies 
(e.g. Giavalisco et al. 2004b). Below z - 3.5, the typical QSO 
colors in the ACS bands move close to the locus of stars and 
low-redshift galaxies. Beyond z — 5.2, the i - z color starts in- 
creasing and infrared bands would be needed to identify QSOs 
efficiently with an "/-dropout" technique. The full set of color 
criteria selects objects in the redshift range (3.5 < z < 5.2), 
while the subset (2-4) is able to select objects with z > 4. To 
avoid contamination from spurious sources, we limited our se- 
lection to Z850 detections with S /N > 5. 

2.2.2. Matching the color selection to X-ray catalogs 

The optical candidates selected with the criteria (1-4) were 
matched with X-ray sources detected by Chandra (A03, G02) 
within an error radius corresponding to the 3 cr X-ray posi- 
tional uncertainty. With this tolerance, the expected number of 
false matches was five and indeed two misidentifications, i.e. 
cases in which a brighter optical source lies closer to the X- 
ray position, were rejected (both in the CDF-S). As already 
shown in paper I, given the flux limits of the Chandra surveys, 
Type-1 QSOs with M145 < -21 should be detectable up to 
z ^ 5.2, up to an optical-to-X-ray flux ratio a ox <: -1.7. We 
compare this result with the estimate of a ox statistics given by 
Steffen et al. (2006, their Table 5). They observe a mean value 
(a M = -1.408 + 0.165) at this optical luminosity, lower than 
our limit. Conversely, any z > 3.5 source in the GOODS region 
detected in the X-rays must harbor an AGN (L x (0.5 - 2 keV) S 
lO^ergs- 1 ). 

The resulting sample consists of 16 candidates, 10 in CDF- 
S and 6 in HDF-N (Table [TJ. With respect to the candidates 
in paper I the following differences are found: in the HDF-N 
(CDF-S) three (one) candidates selected in paper I disappear 
from the selection, while three (one) new candidates enter the 
selection. These differences are due to the improvement both in 
the photometry and in the astrometry of the vl .0 GOODS/ACS 



catalogue with respect to v0.5. In particular, a better matching 
is possible between the optical and the X-ray sources. We dis- 
cuss the implication of these changes in more detail in the next 
section. 

2.2.3. Spectroscopic follow-up 

Spectroscopic information for all 16 candidates has been ob- 
tained by various sources, as listed in Table particularly 
in the framework of the GOODS Team Spectroscopic Survey 
(Vanzella et al. 2005, 2006). 

Thirteen (5 in the HDF-N and 8 in CDF-S) candidates turn 
out to be AGN, with 1 1 QSOs at z > 3 (4 and 7), of which 2 (0 
and 2) are identified as Type II. One object (GDS J033229.29- 
275619.5) in particular turns out to be a QSO with a redshift 
z - 4.759. This candidate was observed in the framework of the 
GOODS/FORS2 spectroscopic follow-up (Vanzella et al. 2006) 
and its reduced spectrum (one-hour exposure with FORS2 at 
ESO - VLT) is shown in Fig.|2ja). In addition to the already 
known QSO GDS J123647.97 +620941. 7 at z = 5.189 (Barger 
et al. 2002), GDS J033229.29-275619.5 brings to two the total 
number of objects in the redshift range 4 < z < 5.2, where our 
selection criteria are expected to be most complete and reliable. 

If we consider the changes in the selected candidates be- 
tween paper I and the present work, we can also draw some 
interesting conclusions. As shown in the previous section, 3 ob- 
jects selected in paper I are not recovered in the present work. 
Two of these objects turn out to be lower-redshift galaxies. 
The third object is a confirmed QSO candidate with a redshift 
z - 2.573. Our analysis shows that its new colors do not satisfy 
our criteria anymore. On the other hand, if we consider objects 
present both in the first and in the present selections, those are 
confirmed QSOs. Finally, 3 new objects are selected with the 
updated version of the ACS-catalogue and all are confirmed 
QSOs. Among them we found an additional candidate (GDS 
J033240.82-275041.4), which is a low-redshift type-2 AGN. 

2.2.4. Checking the completeness of the X-ray 
criterion 

It is worth noting that the z = 4.76 QSO (GDS J033229.29- 
275619.5) is a faint X-ray source that was only selected by 
matching our optical candidates with the supplementary X-ray 
catalogue by G02. This brings up the possibility that our selec- 
tion process might be missing QSOs, in particular at z > 4, if 
their X-ray flux is reduced below detectable levels (e.g. by pe- 
culiar absorption). To assess this potential problem we replaced 
the X-ray criterion (Sect. 12.2.21 with a morphological one and 
applied it to a subset of the color-selected candidates at z > 4. 

The SExtractor software provides a set of parameters quan- 
tifying the radius (in pixels) inside which a given percentage of 
the source flux is concentrated. In the following we fix this per- 
centage at 10%, and we adopt as a diagnostic quantity the re- 
lated parameter (we refer to it as FLUX_RADIUS). SExtractor 
also provides a parameter termed CLASS_STAR, that quan- 
tifies the stellar appearance of the source. For magnitudes in 
the range of interest, point-like sources (CLASSJSTAR ^ 1) 



Table 1. QSO candidates. 
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32 04.93" 


44 31.7 


+ 1.7 


-0.5 


23.668 


1.41 


0.05 


0.20 


< 0.74 


4.10 


1.75 


3.462 


QSO 


Cristiani et al. (2000) 


32 14.44 


44 56.6 


+0.0 


+0.0 


23.094 


2.77 


0.43 


1.42 


0.41 


<4.3 


3.58 


0.738 


Galaxy 


Szokolv et al. (2004) 


32 18.83 


51 35.4 


+0.0 


+0.0 


25.085 


2.26 


0.16 


0.49 


0.74 


14.0 


2.58 


3.660 


Type-2 QSO 


Szokolv et al. (2004) 


32 19.40 c 


47 28.3 


+0.8 


-0.7 


24.676 


1.29 


0.05 


0.42 


< 0.31 


1.17 


1.96 


3.700 


Galaxy 


Vanzella et al. (in orea) 


32 29.29 rf 


56 19.4 


-1.0 


+0.1 


24.984 


> 3 


0.14 


1.72 


0.50 


< 7.6 


1.31 


4.759 


QSO 


Vanzella et al. (2006) 


32 29.84 


51 05.8 


+0.0 


+0.0 


24.642 


2.36 


-0.04 


0.55 


3.06 


31.8 


1.97 


3.700 


Type-2 QSO 


Norman et al. (2002) 


32 39.66 


48 50.6 


+0.0 


+0.0 


24.547 


3.20 


0.21 


0.95 


7.48 


70.6 


2.46 


3.064 


Type-2 QSO 


Szokolv et al. (2004) 


32 40.83 


55 46.7 


+0.0 


+0.0 


25.183 


1.77 


0.09 


0.79 


5.42 


93.4 


2.72 


0.625 


Type-2 AGN 


Szokolv et al. (2004) 


32 41.85 


52 02.5 


+0.0 


+0.0 


22.430 


1.46 


-0.02 


0.48 


16.7 


38.2 


1.25 


3.592 


QSO 


Szokolv et al. (2004) 


32 42.83 


47 02.4 


+0.0 


+0.0 


25.184 


1.43 


0.00 


-0.13 


6.31 


16.4 


1.22 


3.193 


QSO 


Szokolv et al. (2004) 








High-redshift QSO candidates in the HDF-N. 
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8.05 


11.2 


1.21 
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QSO 
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36 42.22 


17 11.6 


+0.2 


+0.2 


24.022 


1.28 


0.23 


0.39 


9.48 


23.2 


1.50 


2.724 


QSO 


Barseretal. (2001) 


36 43.09 


11 08.8 


+0.2 


-0.6 


22.918 


1.11 


0.06 


0.48 


0.75 


< 2.0 


7.13 


3.234 e 




Cowie et al. (2004) 


36 47.96 


09 41.6 


-0.3 


-0.1 


23.761 


4.81 


0.15 


2.07 


2.72 


4.89 


1.27 


5.186 


QSO 


Barser et al. (2002) 


37 03.98 


11 57.8 


-0.7 


-0.8 


25.067 


1.36 


0.05 


0.18 


0.37 


3.49 


1.38 


3.406 


QSO 


Barser et al. (2003) 


37 23.71 


21 13.3 


+0.2 


+0.7 


23.621 


1.21 


0.02 


0.32 


5.07 


11.6 


1.37 


3.524 


QSO 


Cowie et al. (2004) 



" Objects with upper limits both in the soft and hard band have been d etected in other sub-bands. 

* From the supplementary X-ray catalogue by Alexander et al. 1 2003) 

' From the supplementary X-ray catalogue bv^fexander et^y ^^ ^) 

d Detected in the X-rav bvlOiacconi et alJl2002h and not bv lAlexander et alJ<2003h 

e Multiple object: ICohen et alJ feOOOl report z = 0.299. 
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Fig. 1. Left Panel: the B445 - Veeo versus 2775 - z&sa color in HDF-N and CDF-S. Right Panel: the V(,(,o - ins versus 2775 - z^so 
color in HDF-N and CDF-S. The color of point-like candidates are indicated with filled (blue) squares (HDFN) and filled (red) 
triangles (CDFS). The empty circles mark the positions of our 8 final candidates in the combined fields. The solid line shows the 
locus of expected QSO colors at 3 < z < 5.5; we have marked with asterisks the positions corresponding to redshift 3.5, 4 and 
5.2. 



and extended sources (CLASS_STAR a 0) are segregated well, 
while they tend to mix at fainter magnitudes. 

The analysis of the distribution of the FLUX_RADIUS val- 
ues in our sample of QSO candidates shows that objects with 
greater FLUX_RADIUS values turn out to be either lower red- 
shift contaminant galaxies or Type-2 QSOs. This parameter is 
typically lower than 2 for Type-1 QSOs, while Type-2 QSOs 
show FLUX_RADIUS < 2.5. Values greater than 2.5 usually 
correspond to low-redshift galaxies and misidentifications. It is 
interesting to note that, out of the three candidates that were 
selected in paper I and not in the present work (see Sect. 12.2.31 
the two having FLUX_RADIUS > 2.5 turn out to be relatively 
low-redshift galaxies. This criterion can be strengthened if we 
also consider the CLASS.STAR parameter. In this way a com- 
bined morphological selection has been adopted: 

FLUXMADIUS < 2.0 AND CLAS S J> TAR > 0.7 (5) 

To distinguish reliable QSO candidates from stars we compute 
the distance D co /(z) between the point-like candidates and the 
theoretical locus of quasar colors as a function of redshift at 
z > 4, in the (B - V , V - i , i - z) space 



BUZ) = 



ylxLlcol'fizl-colfl 2 

v! 



For B-band dropouts we consider the distance in the (V - 
i , i — z) space: 



Dcol(z) 



V2 



We define the following criterion on the minimum distance 



jynin . 

col 



£>""" < 0.25 mag. 



(6) 



The application of the criteria (5-6) to the ACS catalogues 
produces the sample is listed in Table |3 which consists of 2 
candidates in the northern field and 6 candidates in the southern 
one. Among the candidates in Table |2j we find the two z > 4 
QSOs already listed in Tabled in Fig.^ we plot the colors of 
the candidates listed in Table|5] 

A substantial fraction of these objects were observed dur- 
ing GOODS spectroscopic surveys (Vanzella et al. 2006 for 
the CDF-S, Stern et al. 2006 for HDF-N). The spectra of the 
CDF-S objects are shown in Fig. [2] Two objects in the CDF-S 
(GDS J033240.38-274431.0 and GDS J033218.91-275302.7) 
turn out to be high-redshift galaxies, but they lack any evidence 
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Fig. 2. Spectra of the morphologically+color selected QSO candidates listed in Table|2] 



of AGN activity. One of them (GDS J033218.91-275302.7) 
shows an extremely strong Ly a feature in emission and the 
[NIV] intercombination line, which, together with the absence 
of [NV], is suggestive of an H II galaxy (Fosbury et al. 2003) at 
aredshiftz = 5.554. The spectra of GDS J033201. 28-274432.2 
and GDS J033216.86-275043.9 show no emission lines. The 
quality of the spectra is not sufficient to determine a redshift, 
but enough to exclude their being QSOs in the redshift range 
3.5 < z < 5.2. The proposed object in HDF-N has been ob- 
served in the framework of the GOODS project (D. Stern, pri- 
vate communication). The resulting spectrum is inconclusive, 
but no sign of AGN activity has been revealed. 



On the basis of the morphological selection, no new QSO 
has been found with respect to the X-ray criterion and no 
evidence gathered for an incompleteness in the selection of 



Sect. 12.2.21 The sample listed in Tablets therefore considered 
a reliable estimate of the total QSO space density at z > 3.5. 



3. Estimating the luminosity function of high-z 
QSOs 

In this section we describe the procedure we adopted in order to 
build up the high-z QSO LF from the joint analysis of GOODS 
and SDSS observations. The two samples are complementary 
in terms of the surveyed area and the brightness limit and are 
characterized by different photometric systems, selection cri- 
teria, and spectroscopic completeness. In Fig. [3] we show the 
position of the observed SDSS and GOODS QSOs in the z-m 
space (where m, refers to the i magnitude either in the GOODS 
or SDSS). It is evident from the figure that the two surveys 
cover two non-overlapping regions of this space. 
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High-redshift QSO candidates in the CDF-S. 



Optical 


Z850 


t/ i 
v — l 


l ~ z 


D — V 


v^iass 


Flux 






spectroscopic 


RA 


DEC 


AB 


AB 


AB 


AB 


star 


Radius 


redshift a 




identification 


3'' + m s 


-27° + ' " 


(mag) 


(mag) 


(mag) 


(mag) 


(pixels) 


(pixels) 








32 01.28 


44 32.2 


25.206 


2.19 


0.15 


> 3 


0.75 


1.97 


4.85 






32 16.86 


50 43.9 


25.190 


0.56 


0.25 


2.82 


0.92 


1.39 


3.81 






32 18.91 


53 02.7 


24.561 


2.56 


0.63 


> 3 


0.85 


1.53 


5.554 


EL-Galaxy 


Vanzella et al. (2006) 


32 22.66 


51 37.1 


25.160 


1.49 


0.47 


> 3 


0.98 


1.21 




Star 


Vanzella et al. (in viev.) 


32 29.29 


56 19.4 


24.984 


1.72 


0.14 


> 3 


0.99 


1.31 


4.759 


QSO 


Vanzella et al. (2006) 


32 40.38 


44 31.0 


25.223 


0.53 


0.01 


3.10 


0.78 


1.65 


4.125 


Galaxy 


Vanzella et al. (in Dreo.) 






High-redshift QSO candidates in the HDF-N. 








12'' + m s 


62° + ' " 




















36 47.96 


09 41.6 


23.761 


2.07 


0.15 


4.81 


0.99 


1.27 


5.186 


QSO 


Barger et al. (2002) 


37 11.81 


11 33.5 


25.141 


1.51 


0.11 


> 3 


0.85 


1.60 


4.700 





" Data in italics refer to the photometric estimate of the redshift 
under the assumption of QSO spectrum. 




Redshift 



Fig. 3. Distribution of the observed QSOs (SDSS: triangles, 
GOODS: circles) in the redshift - m, space. 



3.1. Simulated catalogs 

To combine the GOODS and SDSS information in order to es- 
timate the QSO LF, we adopted the approach of La Franca & 
Cristiani (1997, see the appendix for more details on the pro- 
cedure). 

We assume an LF of the form of a double power law: 



#(£l45,£) 



0*(L*) 



(L145 /£*)-<* + (Lus/L*)-P 
or, expressed in magnitudes: 

<D*(M*) 



4>(M 145 ,z) 



100.4(ff+l)(M l4 5-M*) + lQQ.4d3+l)(M li5 -M*)- 



We model a pure luminosity evolution (PLE) 

L*=L* =21) [(l+ z )/3.lf< 

or in magnitude: 

M* = M* =2 l) - 2.5 k z log[(l + 3/3.1] 

and, alternatively, a pure density evolution (PDE) with a power 
law form 

<D* = 0* =23 [(l + z)/3]** 

or, alternatively, with an exponential form 

<D* = $* =2)e W +z) - 3] . 

Given a value for the break luminosity (L* or M*), the 
slopes of the double power law a and /?, the normalization O*, 
and the redshift evolution parameter k z , we are able to calculate 
the expected number of objects from the LF up to a given mag- 
nitude in a given area of the sky. For each object we extract a 
value of absolute M145 and a redshift, according to the starting 
parameters of the LF. We then associate a template spectrum to 
it, randomly chosen from our library (see Appendix A). Using 
the template k-correction and colors at the corresponding red- 
shift, we simulate the apparent magnitudes in the SDSS and 
ACS photometric systems. We also add the effect of photomet- 
ric errors in each band, estimated according to the analysis of 
the photometric uncertainties as a function of the magnitude 
in the DR3QSO and GOODS catalogues. In this way we end 
up with mock SDSS+GOODS catalogues and we then apply to 
them the selection criteria paper I and R02, obtaining selected 
samples of simulated QSOs. 

3.2. Computing the QSO LF 

The advantage of such an approach is related to being able to 
generate and analyze a set of simulated SDSS+GOODS sam- 
ples with the same properties as the real ones in term of se- 
lection criteria, completeness, color, and redshift distribution. 
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We can associate simulated catalogues to each starting set of 
LF parameters. Using multiple realizations for each parame- 
ter set, we are also able to limit the statistical error in those 
simulations, which we can use to estimate the goodness of the 
agreement of a given LF to the observed sample, by directly 
comparing observed and simulated objects. 

Following the approach described in detail by La Franca & 
Cristiani (1997), we divided the z-m, space into bins according 
to the distribution of the observed objects and we defined a 
X 2 statistic by comparing the simulated and observed objects 
in each bin. By using a minimization technique we can then 
compute the best-fit parameters of the LF. 

The functional form of our LF has five free parameters. In 
previous studies they were estimated in various ranges of red- 
shift. Croom et al. (2004) (hereafter Cr4) have found that the 
LF is approximated well by a PLE in the redshift range between 

0. 4 < z < 2.1. We use their estimated parameters at z = 2.1 as 
a starting point for our PDE and PLE models. We try various 
possibilities, differing in the number of parameters fitted/fixed 
and in the type of evolution. As a first attempt we optimized 
only the evolutionary parameter k : (models listed in Table [3] as 

1, 2, 8, 9, 15, 16), then we fit two parameters k z and M+ (mod- 
els 3, 4, 5, 10, 11, 12, 17, 18, 19), and finally we optimized the 
three parameters k z , M* and the faint-end slope a. 

In a recent paper Richards et al. (2005, hereafter Ri5) com- 
piled a sample of quasars using imaging data from the SDSS 
and spectra taken by the 2dF facility at the Anglo-Australian 
Telescope. Their data are in good agreement with the 2QZ re- 
sults of Cr4 at the bright-end, but they require a steeper faint- 
end slope. We also considered models with the faint-end slope 
taken from Ri5, a = 1.45, and the other parameters from Cr4 
(models 2, 4, 9, 11, 16, 18). 

Fan et al. (2001) found a shallower slope at high redshift for 
the bright-end of LF with respect to Cr4 using SDSS commis- 
sioning data. Richards et al. (2006) confirm these findings. At 
Z < 3 they estimate, in agreement with Cr4, a bright-end slope 
P ^ -3. 31, which progressively flattens to become (3 - -2.5 at 
z > 4. To take this possibility into account, we also tested mod- 
els fixing/? = -2.58 and fitting the other parameters (models 7, 
14,21). 

As a final statistical check of the results, we compared the 
simulated and observed QSO distributions with a Kolmogorov- 
Smirnov bidimensional test (see e.g. Fasano & Franceschini 
1987), producing the probability Pks 3 listed in the last column 
of Tabled 

3.3. Results and discussion 

Table [3] summarizes the results for all the models described 
in the previous section, indicating for each of them the fixed 
parameters, the fitted parameters, the corresponding value of 
the reduced^ 2 , and the 2-D Kolmogorov-Smirnov probability 
Pks ■ The best-fitting models are shown in Fig. [5] 

3 It is worth noting that in general the goodness of fit x 1 - wlse an d 
Pks -wise is correlated, but the best-fitting parameters are slightly dif- 
ferent - though compatible within the uncertainties - when minimizing 
the^ 2 or maximizing the P KS 



In general it turns out to be impossible to satisfactorily 
reproduce the observed distribution of objects in the m - z 
plane with PLE models. On the other hand, PDE models show 
qualitatively different results according to the number of fixed 
parameters and the assumed evolution. Models with one free 
parameter (1, 2, 8, 9, 15, 16) are not able to reproduce the 
GOODS and SDSS surveys with the assumed M* = -25.116. 
The agreement becomes satisfactory for PDE models where we 
fit both k z and M* (4, 5, 10, 11, 12). We interpret this result as 
an indication of M+ evolution from z ~ 2.1 to redshifts 3.5 - 4. 
We also notice that models with a steeper faint-end slope agree 
with observations better. This finding is confirmed by models 
with three free parameters: in particular Nr. 13a, which shows 
a relatively high probability for a faint-end slope as steep as 
-1.71. Models with an exponential evolution of the density 
tend to agree with the data better than do models with a power- 
law evolution. 

Models with a shallower slope of the bright-end of the LF 
(as suggested by Fan et al. 2003) are not able to reproduce the 
observed distributions: a bright-end slope as steep as Cr4 is re- 
quired. This result is at variance with Richards et al. (2006), 
we ascribe it to the different estimates of the expected QSO 
colors between them and the present work. As explained in 
the Appendix, we have defined a QSO template library start- 
ing from SDSS spectra in the redshift interval 2.2 < z < 2.25. 
The statistical properties of this library differ from the template 
spectra used by the SDSS team to calibrate their selection cri- 
teria. In particular the mean power-law slope of continuum in 
our library is redder than the one assumed in Fan (1999). As a 
consequence (Fig. IA.4> . we estimate a lower completeness of 
the SDSS color criteria with respect to Richards et al. (2002), 
which is, in turn, the main reason for the discrepancy in the 
estimate of the bright-end slope of the LF. 

On the other hand, the present results also imply that a slope 
steeper than Cr4 and similar to Boyle et al. (2000) is required 
for reproducing the observed data for the faint end of the LF. 

The present observations strengthen the conclusions of pa- 
per I about the space density of high-z quasars and the joint 
evolution of galaxies and QSOs. Simple recipes, which are 
based on the assumption that the QSOs shining strictly follows 
the hierarchical merging structure of DMHs, grossly fail to pre- 
dict the space density of high-z QSOs. Models trying to account 
for the complexity of the baryonic matter behavior with respect 
to DM have to be invoked (see paper I and references therein 
for more details about the feedback effects from AGN and stel- 
lar formation). 

Following the method proposed by Barger et al. (2003), we 
can use our best estimate for the LF (Nr. 13a) to compute the 
QSO contribution to the UV background at those redshifts. We 
calculate the number of ionizing photons per baryon produced 
in the redshift interval by the observed AGN LF, and we convert 
the rest frame 1450 A light to the number density of ionizing 
ph otons us i ng the form of the near-ultraviolet spectrum given 
in BVIadau I i 19951) . We find a ratio between ionizing photons 
and baryons on the order of ~ 0.07 between 3.5 < z < 4.0 
and of ~ 0.02 at 4.0 < z < 5.2. This result is similar to what 
Barger et al. (2003) found for 5 < z < 6.5 QSOs, confirming 
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Fig. 4. Analytical fits to the high-z QSO LF: pure density evolution models with exponential evoution. In all panels the solid line 
shows, as a reference, the model with three free parameters (Nr. 13a in Table[3}. Upper Panel: models with one free parameter; 
dotted and dashed lines respectively refer to models with the Croom et al. (2004) and Richards et al. (2005) faint-end (Nr. 10 
and 11). Lower Panel: models with two free parameters; dot-dashed and long-short dashed lines respectively refer to models 
with the Croom et al. (2004) and Richards et al. (2005) faint-end (Nr. 8 and 9). The position of the filled symbols was obtained 
by multiplying the values of model Nr. 13a by the ratio between the number of observed sources and the number of simulated 
sources. La Franca & Cristiani (1997) demonstrated that this technique is less prone to evolutionary biases with respect to the 
conventional 1 /V max technique. Empty triangles refer to the results of the COMB017 survey (Wolf et al. 2003). 



that the QSO contribution to the UV background is insufficient 
for ionizing the IGM at these redshifts. 



4. Summary 

In this paper we used the SDSS and GOODS databases to build 
an LF in the magnitude interval -21.0 < M145 < -28.5 and in 
the redshift interval 3.5 < z < 5.2. In the first part of the pa- 
per we completed the selection of the QSOs candidates in the 
GOODS fields carried out by paper I. We repeated the origi- 
nal analysis but using an improved version of the optical cat- 
alogues, and discuss the results of the spectroscopic follow-up 
of targets. The main conclusions of paper I on the relatively low 
surface density of faint high redshift QSOs are confirmed and 
strengthened. We also explored the possible presence of "X-ray 
faint" AGNs, not detected in the A03 and G02 catalogues. This 
additional search was based on a morphological analysis car- 
ried out with the ACS optical images, selecting point-like can- 
didates via a direct comparison with the expected QSO colors. 



We followed up the seven most promising candidates spectro- 
scopically finding no additional QSO among them and no evi- 
dence of any significant population of faint X-ray high-z AGN 
with magnitude M145 < -21 in our fields. 

We then combined the QSO sample in the GOODS fields 
with the DR3QSO to study the evolution of the QSO LF. In or- 
der to understand the systematics in the different samples, we 
exploited a method for simulating the combined survey. This 
method is based on a library of QSO template spectra, built up 
starting from observed SDSS QSO spectra in the redshift inter- 
val 2 < z < 2.25. Using this library we were able to predict with 
great reliability the color distributions of observed high-z QSO 
in the SDSS and GOODS photometric systems. Assuming a 
parameterization of the LF and its redshift evolution, we were 
then able to simulate the expected object distribution in red- 
shift and apparent magnitude. From the statistical comparison 
between real and simulated samples, we were able to quantita- 
tively evaluate the agreement of the assumed LF with the data. 
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Fig. 5. Analytical fits to the high-z QSO LF: pure density evolution models with power-law evolution. Symbols have the same 
meaning as in Fig.0] 



We tried different LF models and our results point out that 
PDE models show better agreement with observations than do 
PLE models, even if a luminosity evolution of the magnitude 
of the break, M* U5 , from z = 2.1 to z = 3.5 is required to 
obtain a good match. The most interesting result in terms of 
physical models is related to the faint end slope of the LF. 
Parameterizations assuming a relatively steep faint-end slope 
similar to Ri5 score a higher probability with respect to those 
with a faint-end slope as fiat as in Cr4. If we also tried fitting 
also the faint-end slope, we obtain a best-fit value even higher 
than the Ri5 estimate. A similar conclusion has been reached 
by Hao et al. (2005). Finally, we do not find any evidence of a 
flattening at high-z of the bright-end of the LF with respect to 
the Cr4 value, at variance with the results of Fan et al. (2003). 

We used our best estimate for the LF to compute the QSO 
contribution to the UV background at redshift 3.5 < z < 5.2, 
and concluded that the QSO contribution is insufficient for ion- 
izing the IGM at these redshifts. 

When compared to physically motivated models, the 
present results indicate a suppression of the formation or feed- 
ing of low-mass SMBHs inside DMHs at high redshift. Several 
processes can lead to this phenomenon: photoionization heat- 
ing of the gas by UV background (Haiman, Madau & Loeb 
1999), feedback from star formation (Granato et al. 2004), 
or the interaction between the QSO and the host galaxy. In 



a recent paper, Monaco & Fontanot (2006) studied this last 
possibility in detail, providing hints of a new mechanism for 
triggering galactic winds. In a forthcoming paper (Fontanot 
et al. 2006), this feedback mechanism will be inserted into a 
self-consistent model for galaxy and AGN formation (Monaco, 
Fontanot & Taffoni 2006): the present results for the high-z LF 
will be used to constrain the model and estimate the relevance 
of the joint feedback from supernovae and AGN in shaping the 
redshift evolution of AGNs. 
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" The normalization of the Luminosity Function is in all cases O* = 1.67 x 10~ 6 Mpc~ 3 . 

b Parameters for which no uncertainty is quoted have been fixed in the x 2 minimization. 

c Parameters estimated through x 2 minimization. 

d Parameters estimated through (1 - P KS ) minimization 
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Appendix A: 

A.1. QSO template library 

To build up a QSO template library, we decided to use high- 
quality SDSS QSO spectra at lower redshift. We needed a red- 
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Fig. A.l. Comparison between different estimates of IGM ab- 
sorption. Crosses (green in the electronic version) refer to 
Songaila (2004) ob servations; (ma genta) dashed line refers to 
the prediction of the Madau ( 1995) model (applied to the QSO 
template spectrum of CV90), while the (red) solid line shows 
the applied IGM absorption. 



shift interval for which the SDSS sample had the highest pos- 
sible level of completeness and the continuum of the QSOs 
was sampled in the largest possible wavelength interval from 
the Ly a emission upward. The interval 2.2 < z < 2.25 was 
chosen on this basis. Among the objects in the DR3QSO with 
2.2 < z < 2.25, we selected the higher quality spectra. Our 
final sample consisted of 215 objects. We then estimated the 
rest-frame spectrum of each object, extending the continuum 
the blueward of the Ly a line with the continuum fitting tech- 
nique derived from Natali et al. (1998). We defined several 
"continuum windows" along the spectrum and fit the observed 
fluxes with a power law. We used the resulting parameters to 
estimate the intrinsic QSO spectrum (before IGM absorption) 
blueward of the Ly a line. We used each spectrum in our library 
to compute the theoretical QSOs colors at increasing redshift in 
the SDSS and ACS photometric systems. To simulate the QSO 
spectra at different cosmic epochs, we redshifted our spectra 
up to the redshift of interest, then used a model for the IGM 
to simulate the IGM absorption. To test the IGM model pre- 
dictions at the redshift of interest, we compared them with the 
observations of Songaila (2004). We applied the IGM model 
to the CV90 template spectrum and computed the absorption r 
following the same procedure used by Songaila (2004) on real 
spectra. The result is shown in Fig. IA.1I Crosses refer to the 
observed values of Songaila (2004), while the dashed line rep- 
resents the prediction of the Madau ( 1995) model. It is evident 
from the figure that this model overpredicts the Ly a transmitted 
fraction at any redshift higher than 4, and underpredicts the Lyp 
and Ly y at redshift higher than 3. We used the ratio between 
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observations and predictions to compute an empirical correc- 
tion to the iMadau I i 19951) model at the corresponding wave- 
lengths. The solid line in Fig. IA. II represents the result when 
we applied our modified IGM absorption model to the CV90 
template spectrum. 

A.2. QSOs color prediction 

We then computed the expected mean QSO colors using our 
template library, providing a distribution of colors due to the 
different continuum slopes and strengths of the emission lines 
in observed spectra. 

In Fig. IA. 21 we show our results for the SDSS photometric 
system compared to real colors in DR3QSO. The points refer 
to the color of observed QSOs at different redshifts (see cap- 
tion for more details). In the left panel the dashed line refers 
to the R02 SDSS selection criteria for z > 3.5 QSOs. Our pre- 
dictions of QSO colors are represented with solid lines (dot- 
ted lines are 5% and 95% percentiles of color distribution in 
the template library. The photometric errors are not accounted 
for.). The agreement is quite good. 

We can also predict the redshift evolution of QSO colors. In 
Fig. IA.3l we show the g — r, r-i, i-z color (in the SDSS photo- 
metric system) as a function of the redshift compared with ob- 
served QSO colors. The intrinsic scatter around the mean col- 
ors (i.e. not taking photometric errors into account) is shown 
in Fig. IA.2l (the dotted lines refer to 5% and 95% percentiles 
of the color distribution). Similar results hold for the GOODS 
survey. 



sample, finding a similar result. The result shown in Fig. I A. 41 
(left panel) is different than the analogous plot in Richards et al. 
(2006), who estimated the completeness of their criteria to be 
well above 90% in the whole range of redshift of interest. This 
discrepancy is due to the different QSO templates adopted. In 
particular the SDSS selection criteria are tailored on QSO tem- 
plate spectra whose mean continuum slope f v oc y~ r is "bluer" 
(y = 0.5 + 0.3 Fan, 1999) than the mean slope in our template 
library (y = 0.7 + 0.3). The inferred completeness has a di- 
rect consequence on the shape and evolution of the estimated 
QSO LF. Assuming the Richards et al. (2006) completeness, 
the models with a shallow bright-end (Nr. 7 and Nr. 14) score 
the highest probability, reproducing the result by Richards et 
al. (2006). 



A. 3. Completeness against selection criteria 

It is apparent from Figs. IA.2l and lA.3l that we are able to accu- 
rately predict the expected QSO colors as a function of redshift. 
Our template library can then be used to assess the properties 
of the two independent samples. In particular we focused our 
attention on the completeness of the samples selected using the 
Paper I and R02 criteria. In order to obtain a robust estimate of 
this quantity we applied the selection criteria to our template 
library and analyzed the fraction of QSOs recovered at vari- 
ous redshifts. The result is shown in Figure lA~4l (left panel). 
The SDSS selection criteria shows a lower efficiency in se- 
lecting objects at z ~ 4.4 - 4.5. We tried to investigate this 
issue using the DR3QSO. We considered the redshift distri- 
bution of all objects and compared it with the redshift distri- 
bution of DR3QSO objects satisfying the R02 selection crite- 
ria. The result is shown in Fig. I A. 41 (right panel). We found 
a dearth of z ~ 4.0 - 4.5 objects satisfying the selection cri- 
teria with respect to the total sample. To estimate its signifi- 
cance, we computed the number of object in DR3QSO between 
3.7 < z < 4.0 and 4.5 < z < 5.2 (677) and the number of ob- 
ject in the same redshift interval satysfing R02 criteria (384). 
We then rescaled the DR3QSO redshift distribution by the ra- 
tio between these numbers, and we ended up predicting 210 
QSOs with 4.0 < z < 4.5. Only 164 objects satisfying the R02 
criteria were actually observed, a discrepancy at 3 - <x level. 
We also analyzed the same distribution for the Fan et al. (2003) 



14 



Fontanot et al.: The QSO luminosity function at high-z 




g-r r-i 



Fig. A.2. Color diagrams for the confirmed QSOs in DR3QSO. Stars (red in the electronic version) refer to objects with 3.5 < 
z < 4.0; (green) squares refer to objects with 4.0 < z < 5.2; filled circles refer to objects with z > 5.2; (cyan) dots refer to objects 
with z < 3.5. The solid line shows our prediction of QSO colors in the SDSS photometric system. Dotted lines are the 5% and 
95% percentiles of color distribution in the template library. The photometric errors are not accounted for. Dashed line shows the 
selection criteria of Richards et al. (2002) for z > 3.5 QSOs. 
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Fig. A.4. Left panel: Completeness of selection criteria at various redshifts. The solid line refers to the R02 selection criteria for 
SDSS sources. The dashed line refers to paper I selection criteria in the GOODS fields. Right panel: Redshift distribution of 
DR3QSO sources (solid line). Dashed line represents the redshift distribution of DR3QSO sources satisfying the R02 selection 
criteria. Dotted line represents the redshift distribution of the Fan et al. (2003) sample (multiplied by a factor 10 for graphical 
reasons) taken out of the SDSS commissioning data. 



